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Abstract 

We discuss the mathematical aspects of wave field measurements used 
in traveltime inversion from seismograms. The primary information about 
the medium is assumed to be carried by the wave front set and its per- 
turbation with repsect to a hypothetical background medium is to be es- 
timated. By a convincing heuristics a detection procedure for this pertur- 
bation was proposed based on optimization of wave field correlations. We 
investigate its theoretical foundation in simple mathematical case studies 
using the distribution theoretic definition of oscillatory integrals. 

1 Introduction 

In this paper, we investigate how to carry out tomography directly in terms of 
wavefield measurements. Tomography, in its original form, uses a 'measured' 
wavefront set as input in an inversion procedure which is solely (symplectic) ge- 
ometric in nature, viz. based upon finding bicharacteristics that result through 
a canonical relation in matching the measurement. In 'wave-equation' tomog- 
raphy, one aims at replacing the geometric procedure by a wave-solution proce- 
dure, but keeping the wavefront set of the measurements as the primary source 
of information about the medium. 

Following an embedding procedure to formulate the inverse problem, i.e. in- 
troducing a background medium and incident field and a medium contrast and 
scattered field, we then face the problem of detecting perturbations in the wave- 
front set associated with the scattered (perturbed — incident) field. An intuitive 
choice is based upon correlating the perturbed field with the incident field. We 
will show, by example, that such procedure should be carried out delicately. In 



fact, wc conclude that the perturbation of the wavcfront set can be derived from 
the singular support (of the derivative) of the proposed time correlation. 

The outline of the paper is as follows. We briefly review the microlocal rep- 
resentation of solutions to the scalar wave equation (Section 2). In Section 3 
we introduce the measuring process and its mathematical implementation; we 
describe how the wavefront set of the wavefield propagates through this mea- 
suring process. When we perturb the coefficient function in the wave equation 
(the wave speed) the solution representation will be perturbed. In particular, its 
wavefront set will shift in the measurement-variables cotangent bundle. We for- 
mulate the process of correlating, within the measuring process, the perturbed 
representation with the original representation, and identify how such shift ap- 
pears in the result. It is conjectured that the derivative of the (time) correlation 
at any given measurement position has its singular support precisely at the time 
shift associated with the perturbation of the wavefront set. In Section 4 we give 
examples to illustrate the conjecture. Special attention is paid how to define the 
product of distribution solutions within the correlation process. Finally, in Sec- 
tion 5, we discuss a method of detecting the singular support of the correlation 
in time at any measuring position by means of 'localized' Fourier transforms. 
The procedure defines a criterion to develop wave-equation tomography. 



2 Fourier integral representation of wave 
solutions 

The scalar wave equation for acoustic waves in a constant density medium is 
given by 

(1) Pu = f, 
with 

(2) P = d1 + D c{xfD, 

A{x,D) 

where D = —id^- The equation is considered on an open domain f2 c M" and 
in a time interval ]0, T[. 

We decouple the wave equation into its forward and backward components. 
To this end, we introduce the elliptic operator A{x, D) and its square root 
B{x,D) = ^yA{x, D). Decomposing the field according to 

(3) u± = ^u± ^iB{x, D)-^dtU, 
in combination with the source decomposition 

(4) f± = ±^iB{x,D)-'f, 
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then resuhs in the equivalent system of equations 
(5) [dt±iBix,D)]u±= U. 

Throughout, we assume that c £ C°°{fl). We will construct operators G± with 
distribution kernels G±{x, xq, t, to) that solve the initial value problem equivalent 
to (|) with /± = ±S. 

Let H — H{x,£^,t) — t ± _BP'''"(x,^) denote the Hamiltonian either for the 
forward or backward wave propagation. The Hamilton system of equations that 
generates the Hamiltonian flow or bicharacteristics is given by 

(6) 

dX ~^ dx ' 9A 

Observe that H{x,£„t) = implies r = T-BP™(a:,0- 

Equation can be solved, microlocally, in the form of a Fourier integral rep- 
resentation. The phase of the associated Fourier integral operator follows from 
the canonical relations 

C± = {(a;(xo,eo,±i),t,e(a;o,eo,±t),TSP""(xo,eo) ; xo,-^o}- 

. ■' 

T 

Let 

{xi,XQ,ij,T) with / U J = {1, . . . ,«} 

e 

denote coordinates on C± . A function S will locally describe C+ according to 

(7) ^ ^ 

= — S , ^0 = -- — S , 

OXi OXq 

and generates the non-degenerate phase function 

(8) (t)+{x,XQ,t,£,j,T) = S{xi,xq,^j,t) + {^j\xj) +Tt. 

In our notation, we will suppress the dependence on xq and collect £,JtT in the 
phase variables 9. The canonical relation can then be written as 

C+ = {{{x,t,d^(j)+,dt(j)+); {xo,-dx„(f>+)) \ de(t)+ = 0}. 

We synthesize the canonical relation — C+ U C_ with associated (non- 
degenerate) phase function = 0_ if r > 0, = (/)+ if r < 0. In accordance 
with (|) we obtain 

(9) G{x,xo,t) = ^i[G+{x,xo,t)-G-{x,xo,t)]B{xo,D,„)-\ 
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With this fundamental solution, the solution of (||) and its dependence on the 
initial conditions can then be written in the form of a Fourier integral operator 
(FIO) with amplitude a = a{xi,xo,S,j,T). In fact, a is a section of the tensor 
product Mc^ 0f2^/^(C0) of the Keller-Maslov line bundle and the half-densities 
on Cm- 

The kernel of the FIO admits an oscillatory integral (01) representation. In 
the remainder of this paper we consider such OIs to represent 'the wavefleld'. 
Perturbation of this wavefield are induced by perturbation of the coefficient 
function c(x). 



3 Detection of singularities of the wave field 



As described above, each component of the wave field as well as the perturbed 
wave field can be represented by an 01, 



(10) 



u{x,t) = / e''^(^'*'^'a(a;,t,6l)d6l, 



where (p is a non-degenerate phase function and a a symbol (| Hor90 |, Sect. 7.8); 
note that the wave front set satisfies the inclusion (| Hor9Ct| , Thm. 8.1.9) 



(11) 



WF(u) C {(x, t; a,0(x, t, 9), dt^ix, t, 9)) \ dg^ix, t, 9) = 0}. 



3.1 Measurements as restrictions to submanifolds 

Measurements are recordings of the wave field u in stations at certain points 
x in the acquisition manifold over some time interval {to,ti); mathematically, 
this corresponds to the restriction of the distribution u to the one-dimensional 
submanifolds Sx = {x} x R followed by further restriction of the resulting one- 
dimensional distribution Ux of time to the open interval (to, ti). 

While the second of those restrictions is always possible and straightforward, 
the first can be carried out as continuous map only on distributions satisfying 
the following condition ( |Hor90| , Thm. 8.2.4 and Cor. 8.2.7) 



(12) WF(w) n {{x, i; 0) 1 1 e M, M"} = 0. 

Note that by ( pT| ) this condition is satisfied if and only if dt4>{x, t,9) ^ when- 
ever dg(j){x, t, 9) = 0. If it holds, the restriction can be defined as the puUback 
of u under the embedding map '■ Sx ^ R"+^ and by ( pi] ) we have the 
wave front set relation 

(13) WF(u,) C {{t- r) I 3^ (x, t- r) S WF(m)} 

C {{t-dt(P{x,t,9)) I de<P{x,t,9)^Q}. 
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Let ^ be a phase function, b a symbol, both with the same domains and supports 
as 4>, a, and v be the oscillatory integral defined by them; assume that v also 
satisfies (|l^) and set Vx = i^v. In case we are interested only in a certain time 
window of measurement we may use further cut-offs and achieve that and 
Vx are compactly supported. 



3.2 The correlation function 

For a e R denote by Ta the translation by a on R. If the distributional product 
Wx^t = 'u-x • TiVx can be defined and yields an integrable distribution ([ Hor66 |, 
Sect. 4.5) we define the value of the correlation Junction at t by 



(14) 



c[Ux,Vx\{t) = {Ux 



TtVx 



I) 



{Wx.t, 1) 



The correlation is bilinear in [., .]. Whenever there is no ambiguity about the 
distributions u and v and the point x under consideration we will denote the 
correlation briefly by c{t). 

Whenever u and v represent the unperturbed and perturbed solution, then 
typically dt<t> = dti/j (the frequencies coincide) and therefore for certain values 
of t we expect the cotangent components of the wave front sets of Ux and TfVx 
to be identical on the overlap of singular supports. That means that, unless 
both cotangent parts are only half rays on the same side of 0, Hormander's 
condition ([Hor9C|, Thm. 8.2.1) for defining the product does not apply. But 
within the hierarchy of distributional products described by Oberguggenberger 
(|Obc92], Ch. II) this condition, 'WF favorable', appears only as one out of a 
variety of consistent possibilities to give a distributional meaning to the product 
under consideration. We apply some of these to the analysis of the correlation 
function in some examples below to explore and illustrate whether and how the 
correlation, after restriction, can provide information about shifts in wave front 
set from u to v. It will become clear that the customary criterion of searching 
for the 'stationary point' of the correlation (Dahlen, Hung and Nolet |DIINOO|, 
Zhao, Jordan and Chapman [ZJCOC| and Luo and Schuster |LS91]) for detecting 
the shift in wave front sets is generally incorrect. 

Here, we would like to point out that the appropriate mathematical framework 
to deal with the multiplication (and also the restrictability) in a unifor m and 



system atic manner is Colombeau's theory of generalized functions (cf. [ Col85 



Obe92 |). Such framework will enable us to c ope wi th the integrability question 
(forming the correlation) at the same time (|Hor9£]). 



Practically, we will have to consider regularizations or approximations to the 
formal expression c{t) = {wx,t, 1) of the correlation either to give a meaning to 
the product or to make the integration (i.e., distributional action on 1) well- 
defined. This amounts to the attempt of defining c(t) as the pointwise (in t) 
limit of sequences 



c„(i) = «,,l) 
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as n — > oo where w"^ is a suitable regularization or approximation of Wx^t- 
3.3 The shift of singular supports 

We compare the singular supports, or rather the wave front sets, of Ux and Vx- 
Their offset expresses the amount of time shift of the wave fronts (or rather 
singularities) at location x by the perturbation. 

First we observe that under a natural time evolution condition on the phase 
function a restrictable 01 is representable as an 01 in one dimension. 

Lemma 1. If u = / a(., 6l)e''^(-^) e P'(M"+i) satisfies condition at x 
and dt(f>{x,t,9) / for all t and 6' 7^ such that {x,t,6) £ supp(a) then 
the restriction Ux to Sx is the 01 on M (i.e., in the time variable) where x is 
considered as a parameter in the phase and amplitude. Therefore 



Proof. By assumption (px{t, 9) — (p{x, t, 9) defines a phase function on M x M^. 
We have Ux = tj(u) and t* is continuous on the subspace of restrictable distri- 
butions. Therefore we may use any standard 01 regularization u — lim£_,o 
and obtain Ux — lim^^o '-xl'^^e)- Since the latter is an 01 regularization in one di- 
mension with phase function (f>x and symbol a{x, ., .) the assertion is proved. □ 

Note that the usual stationary phase argument applied to this one-dimensional 
01 gives the same upper bound for the wave front set as established above in 
(|3|). Assuming that the perturbed solution v is given as an 01 with phase 
function ip and amplitude b we can compare the wave front sets of their mea- 
surements at X (restrictions to Sx)- 

As pointed out above, the perturbation will affect the phase function only in 
its X- and ^-gradient, i.e., we may assume that dt4> = dt^j. If {to,To) G WF{ux) 
then To = dt(l){x, to,9o) for some with dg(j){x, to, 9o) — 0; similarly if (ti, n) G 
WF{vx) then n = dt(j){x,ti,9i) for some 9i with dg<j){x,ti,9i) — 0. In any 
microlocal representation of the solution to the wave equation, in the absence 
of attenuation, the phase contains t only linearly, say, in the form tp for some 
conjugate (frequency) variable p. 

As was shown in Section 1, typical phase functions are of the special form 
(j){x, t, 77, p) = (poix, rj, p) — tp and '4>{x, t, 77, p) = ipoix, rj, p) — tp. In this case the 
stationary phase conditions (see eq. (|ll|)) in the wave front sets read 

(16) io = 9p0o(2^,?7o,Po), i9i,0o(2;,'7o,Po) = 

(17) ti = 9p?/'o(2^,?7i:Pi)i dn'4'o[x,rii,pi) = Q 

and the respective t-derivatives of the phases yield cotangent components tq = 
— Po and Ti = — pi. By the (positive) homogeneity of (p and w.r.t. (77, p) 



(15) 
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their first order derivatives w.r.t. those variables are (positively) homogeneous 
of degree 0. 

Hence, if we are detecting time-like forward (resp. backward) cotangent direc- 
tions, i.e., p > (resp. p < 0), we may rescale the arguments in the phase and 
obtain the time shifts 

(18) ti-to = dptpoix,r][,±l) - dp(t)o{x,r]'Q,±l) 

for the corresponding (slowness co- vector) projections 770, satisfying the con- 
ditions 

(19) d^Mx,rio,±l}^0, dr,Mx:V'i,±'^)^0. 



3.4 Correlation optimization 



In [LS91] a traveltime inversion method is described that uses optimal fitting 
of traveltimes from synthetic scismograms according to wave equation solutions 
of velocity model perturbations. The fitting criterion is based upon a crosscor- 
relation function of the observed {v) and the synthetic (m) seismic data. This 



crosscorrelation of |LS91| corresponds to the correlation function defined in ( |14[ ) 
above. 

We give a brief schematic description of this interesting fitting strategy and test 
its theoretical validity in three simple examples below. Assume that v represents 
the observed (or perturbed) wave field and u = ^[7] is the solution of a velocity 
model which is parametrized by the variable velocity j{x). We assume that 7 
is a real-valued smooth function. Therefore the correlation function is actually 
dependent on time t and the velocity 7 which we indicate in the notation 



where (.) denotes the scalar and [.] the functional argument of c. An intuitive 
expectation would then be that at the exact traveltime shift induced by the 
perturbation, we find optimum match (overlap) of the corresponding seismo- 
grams and therefore the crosscorrelation should be maximal. Leaving possible 
maxima at time interval boundaries aside, we search for a (7, t) relation that 
gives stationarity of the crosscorrelation, i.e., 

(20) F{tM^dtc{tM^O. 

Naively speaking we can consider this to be an implicit definition of a functional 
relationship between 7 and t. (Observe that 7 is an infinite-dimensional variable 
and therefore more attention is to be paid to the exact meaning of applying an 
'implicit function theorem' below.) Under the condition that dtF = d^c ^ we 
would therefore try to solve equation (^o|) locally for i as a function of 7 and 
find a quasi-explicit representation by 
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4 Case studies 



4.1 Two propagating delta waves 

Consider u = 5o(x — s) and v — So{x — 7s), two Dirac deltas travelling along 
the lines x = s and x = js respectively. (These are distributional puUbacks 
of 60 G V'(R), the Dirac measure located at 0, via the maps {x,s) t-* x — s 
and (a;,s) 1^ x — 7s.) Assume that x > 0; the opposite sign case is completely 
symmetric. We clearly have Ux = Sx and T^Vx = -^S^-t, and therefore 

X 

singsupp(ua;) = {x}, singsupp(u2;) = {-} 

7 

yielding a singularity shift of <i — io = —x(l — I/7). 



Observe that Ux and T^Vx have disjoint singular supports unless t — —x{l — 
1/7) in which case their product would require to multiply 6x with itself. This 
cannot be done consistently within the hierarchy of distributional products (cf. 
[ 3be92| ) and calls for a systematic treatment in the framework of algebras of 



generalized functions. However, here we touch upon those aspects only in terms 
of regularizations. 

Choose a rapidly decaying smooth function p on M such that / p = 1, in other 
words /9 is a mollifier, and set pds) — p{s/e)/£. Denote by and the 
convolutions of Ux and Vx with p^. Then we have 

1- 



= Peis - x), T*V%{s) = -peis + t - X h) 

1 



and upon integration of u':j.{s)T^v^{s) w.r.t. s with a change of the variable 
y = (s — x) I £ we obtain for the regularized correlation function 



(21) cM = - / p{s)p[s + ^±^(LJZ2:)) ds. 



76 

If we let e — > we observe that c^(t) pointwisc for t ^ — a;(l — I/7) 
and |c£(f)| — > cxD. Hence, in an approximative sense, the singular support of the 
correlation c(i) contains the time shift information. To be more precise, it is 
not difficult to show that in the sense of distributions 

(22) Ce \^8_x(\-il^). 

For this, we just note that for arbitrary Lp e I?(R) one may change the variable 
in J (p{t)cg{t) dt to r = {t — t)/e and use the fact that J f * g = J f ■ J g for 
rapidly decreasing functions / and g. 

In particular, this shows that here the correlation is stable under changes within 
the chosen class of regularizations since the limit does not depend on p. 
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Curiously enough, the regularization approach also gives the correct answer 
when using the procedure of |LS91|. Define the short-hand notation kg{x, t, 7) = 
{t — t)/e and consider 

= — 2 I Pis)p'{s + ke{x,t,j))ds 
IS J 

and set Fs{x,t,'y) — je^c'^{t). We see that 



dteFe{x,t,^) ^ J p{s)p"{s + ke{x,t,j))ds, 

which is proportional to ||p'||l2 at i = t and stays nonzero for all t and 7 close 
enough. This in particular true at 7 = 1 in which case i = 0. Therefore, in 
a (possibly smaller) neighborhood of these values for t and 7 we can solve the 
implicit equation F^{x, i, 7) = for t — t{x, 7) and find locally 



d t(x d^F,{x,t,-i) _ p[s)p"{s + ke[x,t,-i))ds _ x 

dtF,{x,t,^)~ ^ J p{s)p"{s + k,{x,t,j))ds 7'' 

We find from this by integration over 7 (close to 1) that 

t{x,j) = X — ~a;(l — 1/7) 

7 

which is the correct shift of the singular support. 



4.2 A delta wave interacting with a shock 



We set u = 6o{x — t) and v = H{x — -ft) (where H is the Heaviside function) 
yielding exactly the same configuration of wave front sets as in the previous case. 
In this case, restricting our attention again to x > 0, — 5x and T^Vx{s) = 
H{x — 7(5 + t)). The only critical product appears if t = t = —x{l — I/7): At 
this point we have to deal with 6x{s) ■ H{x — s) which exists as a so-called 'strict 
product (7.4)' in the notion of |Obc92|, Ch. II, assigning the value ^Sx to it. 
For t < t we obtain Ux ■ T^Vx = Sx and for t > t we have Ux ■ T^Vx = because 
the Heaviside contribution is constant or 1 in those regions. In summary 



Wx,t = Ux ■ TTVx = < 



5x 


iit<t 








v 


\it>t 



If we interpret {wx,t, 1) via the Fourier transform of Wx,t as w^ti^^) 
we obtain c{t) as the measurable function 



c{t) 



c{t) then 
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As in the previous case, we observe that it is exactly the singular support of c 
- here, the point t = — x(l — I/7) - that reveals the information of the correct 
shift. 

Observe, however, that the travel time t is in fact the only point where the 
(distributional) derivative c'(t) = St does not vanish. The previous evaluation 
based upon the implicit function theorem hence does not apply. 



4.3 Wave equations with different medium constants 

We now return to the wave equation (Section 2), assume constant coefficients, 
and invoke an exact solution representation rather than an asymptotic one. We 
consider propagation in one spatial dimension. 

Let X £ C°°(M) be real- valued, x{~0 — xiOi X = ^ neighborhood of and 
X = 1 for 1^1 > 1; 7 a constant > 0. 

In the sense of OIs 

u{x,t) = 

vix,t) = /ei(*l«l-C/7)X(eZ7)^^ 

(it respectively v are the complex conjugates of 27ri times the subtrahends in 01 
representation of the fundamental solutions for the d'Alembert operator with 
wavespeed equal to 1 and 7, respectively.) 

The general WF-bounds according to (|l^) give 

WF{u) C {{t,±t,~t\i\)\teR,±^>0} 
WF(z;) C {(7t,±t,-e/7,|ei) |te]R,±e>0}. 

Observe that half rays in cotangent components are minimal closed cones in 
\ 0. We will show that, in fact, the inclusion should be replaced by equality. 

For symmetry reasons, we give detailed arguments in quadrant x > 0, t > 
only. Since {df — d'^)u — and [df — j'^d'^)v — 0, the theorem on the 



propagation of singularities |Hor90|, Thm. 8.3.3, applies; in particular, if {t,t) G 
singsupp(u) (resp. (7^, i) G singsupp(i;)) then the whole line through this point 
with directional vector (1,1) (resp. (7, 1)) is in the singular support (with the 
same perpendicular cotangent component in the wave front set attached to it). 
Therefore, to prove equality in the above WF inclusion relations, it suffices to 
show that u (resp. v) is not smooth near (0, 0) 

Assuming the contrary, would imply that the function x t—^ dtu{x, 0) is smooth; 
but it is also equal to the Fourier transform of 1%, which cannot be smooth since 
J-'X = T{x — 1 + 1) = J'ix ~ 1) + 27r(5 where the first term is a smooth function 
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(of rapid decay) since x — 1 is smooth and of compact support. The argument 
for dtv{x, 0) is the same. We conclude that 

(23) WF(it) = {(t,±t,-e, lei) |t eIR,±e>0} 

(24) WFiv) = {(7t,±t,-e/7,iei) UeK,±e>0}. 

It foUows immediately that both u and v are restrictable to Sx and and Vx 
are represented as the one-dimensional OIs where x appears as parameter in the 
phase only. Assuming a; > 0, we clearly have 

(25) WF(u,) = {(±a:,OI^>0} 

(26) WF{vx) = {(±a;/7,0 U>0}. 

But then the time shift is given by 

(27) ti - to - T^il - 1/7), 
as expected from physical intuition. 

In the remainder of this section we analyze the correlation in detail, and inves- 
tigate how the time shift appears. In the correlation we have to multiply the 
distributions 

u.(.)=|e*W«l-«)^d^ 

and 

J 7l?l 

which have wave front sets 

WF{ux) = {-a;,a;}xM+ 



WF{TtVx) = { t, - - t} X M+. 

7 7 

Hence, whenever t ^ T2;(l T 1/7)1 the distributions have disjoint singular sup- 
ports and in case t = ^x{l =F I/7) the cotangent vectors in their wave front 
set cannot add up to 0. We conclude that for all t the wave front sets are in 
favorable position and the product w^.t = Ux ■ TtVx G I?'(M) can be defined in 



the sense of [ Hor9C | , Thm. 8.2.10. The following lemma states that we are even 



allowed to use the naive product of the 01 expressions. 
Lemma 2. Wx^t is (essentially) an 01 given by 

(28) wxt{s)^ II e'(^(l«l-l''l)-*l''l--(«-''/^)) d^dr^ 

JJ iiQim 
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and t I— > Wx,t is weakly continuous R 2?'(R). Wc introduce the following 
notation: 

(l>t{s-A,v) = s{\£.\-\v\)-'t\v\-x{^-v/i) 
i\Q\v\ 

for the phase function and the amplitude. 



Proof. For the just ification of we use the construction of the distributional 
product in Hor9C| , Thm. 8.2.10 via the pullback of the tensor product on 



2 



under the map l{s) = (s, s) which embeds K as the diagonal into M . In doing 
so the original OIs may be approximated by smooth regularizations (e.g., am- 
plitude cut-offs in the integrands) the tensor products thereof being pulled back 
simply as smooth functions (meaning restriction to (s, s) in this case). 

It is easily seen then that the smooth functions obtained thereby converge weakly 
(as 01 regularizations) to the 01 given in (|28|). By continuity of the pullback 
(under the given wave front set conditions) this limit equals the pullback of 
the tensor product of the corresponding limits and therefore, in turn, is the 
distributional product Wx,t — ■ Ttv^. 

Note that a is smooth in (^, rf) (due to the cut-off x) and homogeneous of degree 
—2 outside the set {|^| > l,\ri\ >7} and is therefore a symbol of order —2. The 
function 0t is smooth on supp(a) and homogeneous of degree 1 in (^,?7). If 
\t\ ^ |a;(l± 1/7)1 then the gradient d(s^^^ri)4>t 7^ (0, 0, 0) for all (s, ^, 77) and hence 
</>t is a phase function. 

In case |t| = \x{l ± 1/7)] the gradient vanishes exactly along one half-ray com- 
ponent of the set {(1^,77) | |C| — \ri\} (e.g., along ^ = ?7>0ift = —x{l — I/7)). 
Although it is no longer a phase function in the strict sense, the distribution 
Wx^t is then defined as the sum of a classical integral, an 01, and a Fourier 
transform of an L^-function. We discuss this for the case t — —x{l — I/7) in 
detail, the other cases are completely analogous. 

Let /z(^, 77) be a smooth function that is equal to 1 near ^ = 77 > 1, has support 
in > 0,77 > 0}, and satisfies < /i < 1. Let 2^(^,77) be smooth with compact 
support and 7^(^, 77) = 1 when + 77^ < 1. 



12 



(l-v)(l-^) 




We can split the integral defining Wx.t into three terms according to 1 = i>{l — 
n) + {1 — — fj.) + ^. The first integral, then, is a classical one defining a 
smooth function, the second is an 01 since the gradient of </>( does not vanish 
on the support of the integrand. In the third integral, we have 



(pt{s; ^, ?7) = (s - x){^ -v) = -{{x 



(insert t — ~x{l — I/7) and use the fact that |^| = C > ^^^d jryj = > on 
the support of the integrand) and hence the last term is equal to 



which we interpret via the Fourier transform of the L^-function fm on as 
s '—^ J-'{fia){x — s, s — x) in the sense of locally integrable functions — hence it 
is distribution on M. 

The weak continuity w.r.t. t follows from the smooth dependence of the phase 
function in the 01 representation (cf. |Dui96|, before Thm. 2.2.2) and the con- 
tinuity of the Fourier transform on L^. □ 



Remark 3. From the last part of the proof it follows that t (-^ Wx^t is weakly 
smooth on M \ {±a;(l ± I/7)}. 

In order to define the correlation function, we need to check whether the action 
of Wx^t on 1 is well defined. We will do so by showing that Wx^t is tempered 
with Fourier transform w^t being in fact a continuous function. This function 
can be evaluated at yielding the interpretation {wx^t, 1) = WxTtiO)- 

We use an 01 regularization of Wx,t via the symmetric cut-off function p(S^, rj) = 
PoiOPoiv) where pQ G I?(R) with po{r) — when \r\ > 1, po{r) = 1 when 
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\r\ < 1/2 and < po < 1. Writing p.i^T]) = p{^/j,v/j) = 1,2,...) we 
obtain supp(pj) C [— and pj 1 uniformly over compact subsets of as 
j ^ oo. Hence 



w. 



,,t = D' - lim / / e''^*(-«''')aj(^,r?) d^dr] = V' - lim wi t 



where 

Since supp(aj) C [— j, j]^ is compact s 10^.^(3) is smooth and by differenti- 
ating inside the integral we see that for all I G No ij^Ywi t(*) bounded by 
some constant (depending on I and aj). Hence (w;^ J^gN is a sequence in the 
space 4S'(M) of tempered distributions. Therefore, to prove that Wx.t is in tS'(M), 
it sufSces to show that (w^ f)ieN converges weakly in 5'(M), i.e., for all rapidly 
decaying smooth functions ip £ <S(R) the sequence (w^ j, ip) is convergent. 

We have 

= J e-'(*l''l+^(«-''/T»a,-(Or/)ye^«(l€l-l''l)v?(s)dsrf(C,??) 



Here, the integrand tends pointwise to e~^^*^'^^~^^'^^~^^'^^'>a{^,ri)(p{\r]\ — |^|) as 
j —> 00 and is dominated by \a{^,r])(p{\ri\ — |5|)|. It remains to show that 
(?,^) o(0?7)^(|r7| — 1^1) is in L^(M^); then an application of Lebesgue's domi- 
nated convergence theorem will provide us with existence of an explicit integral 
expression for the Hmit {wx,t, v)- 

Since ^ € '5(M^), using the explicit structure of a, we have for any fc e N a 
bound of the form 

Hi,v)Q{H - m < Ck{i + + \v\)-\i +H- mr' e R'. 

While integrating the right-hand side of this inequality over R^, wc split the 
integration into four parts according to the sign combinations of ^ and r]. By 
symmetry, this boils down to estimating only the two kinds of integrals 

CXD 00 00 00 

d^drj f f d^drj 



(l + 0(l + r?)(l + h-C|)'=' ^ J J {l+Oi^ + v)a + V + 0'' 
00 00 
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In /+ wc only have to note that (1 + C + v)~'' < (1 + + f?)"''^^ which 

together with the remaining factors gives a finite integral as soon as fc > 0. In 
/_ we change variables to = 77 — ^, fi — ?]to obtain 

00 00 

/-= / — ^ / . .d.. 



max(0,z^) 



In the inner integral we use 1 + ji — iy — (l + /i)(l — + > (l + /i)/(l + |i^| 
yielding an upper bound (I + I^^D/q + and hence 



d/i f dv 



which is finite if fc > 2. This proves the assertion that 77) ^ a(^, rf)^{\r}\ ~~ |^|) 
is indeed in L^(M^) and establishes the following result. 

Proposition 4. Wx,t G 5'(R) and for any ip G iS(M) 

(29) {w^.uv) = lim K.*,.^) = / e-'(*l''l+-(«-''/^»a(e, r7)^(|77| - \(.\) d{^,r^). 

We are now in a position to determine the Fourier transform of Wx.t explicitly. 

Proposition 5. w^t is the continuous function on M given by (the classical 
integral) 

(30) .-,M=4„*/e-<*-ICI.eo.,i,|£|-.„o,«,K| 

{I«l>'-} 

Proof. Let Lp G iS(M) then (p{s) = 2inp{—s) and from ( |29| ) we obtain 

= 2^1 e-*(*l''l+-«-''/^))a(C, 77)9,(1^1 - \v\)d{tv) 




= 2tt J e-"i(^ J e-'^-''^-^'^/^^ai^,r])ip{\^\+r])dr, 

— 00 

+ I e-'<'^--^/'^'>a{tv)^m~T^)dT^)dt 



where in the last line we have made use of the symmetry properties of a{^^rj). 
Changing coordinates in the inner integrals to r = |^| ± 77 and again by the 
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symmetry of a(^, 77) this reads 

2tt /"e-'(^«+*l«l) / e"''a(C,|<e|-r)95(r)(e'('^-l€l)^/'^-e-'(''"l«l)="/'') drd^. 



-°° 2cos(^(|5|-r)) 

Finally, since e 5 and \a{i, \i\ - r)\ < p{r){l + for some polynomial in 

r, we may interchange the order of integration and arrive at 

{^uV)^ J ifir) ■ 4ne^"- J e-'^-^+*\^\^ cosi-^ ~ r))a{^,\^\ - r) dr. 

m>r} 

Since tp was arbitrary and the above upper bound for a(^, |^| — r) shows that 
the inner integrand is in w.r.t. ^, the proposition is proved. □ 

From ( |30| ) we immediately obtain the correlation by setting c{t) — w^t{0), in 
the form 

00 

(31) c(t) = 47r y" e-'("«+*l«l)cos (^^^ a(e,0'^^ 



00 

t(/ 



00 



-i(xC(l~l/7)+t|?|)^^^iiMll 

„-i(i.5(i+i/7)+t|ci):^iiMil 



-00 



This shows that t <—>■ c{t) is continuous and can be represented as the difference of 
two (classically convergent) OIs with symbols of order —2, and hence c £ L^(R). 
Note that the (distributional) derivative c'{t) can be obtained by differentiating 
w.r.t. t inside the 01 raising the order of the symbol by one. Therefore c' will 
not be continuous on the whole line. 

Finally, we observe that again the information about the singularity shift is 
revealed by the singular support of c{t). By the stationarity condition on the 
phase functions, we find 

WF(c) C {±x(l + 1/7), ±a;(l - I/7)} x R+, 

where ±x(l — I/7) represent the true shifts from ±x/j to ±x whereas ±a:(l + 
1/7) are the distances from T^/j to ±x. It is easily seen that c{t) cannot be 
smooth at the points t = ±a;(l ± I/7), e.g., by noting that each time derivative 
brings down a new factor of |^| in each integrand, and at the t values in question 
one of the phase functions vanishes identically along a half-line in ^. Hence, we 
have in fact the exact information 

(32) singsupp(c) = {±x(l + I/7), ±x(l - I/7)}, 

which also fits nicely with remark |^ on the weak smoothness of Wx,t- 
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5 Microlocalization of the correlation 



From the case studies, we conjecture that the singular support of the correlation 
of two wave fields reveals the relative shift in wave front sets between them. As 
we pointed out, in general, the critical point set of the correlation need not be 
compatible with this shift. Here, we propose an alternative approach to extract 
the shift from the correlation, viz., by detecting its singular support. We design 
a pseudodifferential operator that enables this detection. Our approach can be 
applied invariably to any derivative of the correlation also. 

In the generic case, the correlation c £ S' with Fourier transform c. Let </) be 
the Gaussian in one dimension, define 

1 f s 

(33) i'rA'^) = -4> 

r 

Introduce 

(34) W^^,c{^) =i^t*c{T) for T = ±l, 

a continuous wavelet transform that can be written as the action of a pseudod- 
ifferential operator V'r.tl-Dr) (in Op S~°°) on c. The growth properties reveal 
the wave front set at t in the direction ±1. In fact, (i, ±1) ^ WFc if for any 

N €N, 

(35) |Wv,,,.c(±l)| < CArr^ for re]0,l] 



(see [lFol89[ ). Effectively, this leads to a scanning procedure over t: whenever 
the condition is not satisfied, t G singsupp(c). In particular, this applies if 
|W^^,c(±l)| w r*^ for some fixed M. 

If c would allow an 01 representation, as is the case in the examples of Section 4, 
we could apply a stationary phase argument instead, as in (hj 



6 Discussion 



Starting from the microlocal representation, we analyzed the measurement pro- 
cess of wave fields. Such process can be described by a restriction operator. 
We then adressed the issue of how the detection of wave front sets propagates 
through the measurement process. Then wc focused on the detection of (base) 
shifts in wave front sets due to perturbation of the wave field within the mea- 
surement. We introduced the distributional cross-correlation as a tool for this 
purpose, and analyzed its properties. In a series of case studies, we investigated 
in what way the cross-correlation reveals the shifts. In the first case the corre- 
lation was a measure, in the second case it was a bounded measurable function, 
and in the third case it was a continuous function. It was conjectured that the 
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time shift coincides with the singular support of the correlation. We proposed a 
procedure (a pseudodifferential operator) to detect the shift based on microlo- 
calization. Such procedure would comprise the foundation for wave-equation 
tomography. 

Acknowledgement: We thank Jerome Le Rousseau for valuable mathemat- 
ical remarks and improvements of the text. 

References 

[Col85] J. F. Colombcau. Elementary Introduction to New Generalized Func- 
tions. North- Holland Mathematics Studies 113. Elsevier Science Pub- 
lishers, 1985. 

[DHNOO] F.A. Dahlen, S.-H. Hung, and G. Nolet. Frechet kernels for finite- 
frequency traveltimes-I. Theory. Geophys. J. Int., 141:157-174, 2000. 

J. J. Duistermaat. Fourier Integral Operators. Birkhauser, Boston, 
1996. 

G.B. FoUand. Harmonic analysis in phase space. Princeton University 

Press, Princeton, 1989. 

J. Horvath. Topological vector spaces and distributions. Addison- 
Wesley, Reading, MA, 1966. 

L. Hormander. The Analysis of Linear Partial Differential Operators, 
volume I. Springer- Verlag, second edition, 1990. 

G. Hermann. Integration and microlocal analysis in Colombeau alge- 
bras. J. Math. Anal. Appl, 239:332-348, 1999. 

Y. Luo and G. T. Schuster. Wave-equation travel time inversion. 

Geophysics, 56:645-653, 1991. 

M. Obcrguggcnbcrgcr. Multiplication of Distributions and Applica- 
tions to Partial Differential Equations. Pitman Research Notes in 
Mathematics 259. Longman Scientific & Technical, 1992. 

L. Zhao, T.H. Jordan, and C.H. Chapman. Three-dimensional Frechet 
differential kernels for seismic delay times. Geophys. J. Int., 141:558- 
576, 2000. 



[Dui96] 

[Fol89] 

[Hor66] 

[H6r90] 

[H6r99] 

[LS91] 

[Obe92] 

[ZJCOO] 



18 



